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We observe a new D s meson with mass (2856.6±1.5 s tat. ±5.0 SW st.) MeV/c 2 and width (48±7 s t a t. ± 
Wsyst.) MeV/c 2 decaying into D°K + and D + K° s . In the same mass distributions we also observe 
a broad structure with mass (2688 ± A sta t. ± Ssyst.) MeV/c 2 and width (112 ± 7 sta t. ± 36 syst .) 
MeV/c 2 . To obtain this result we use 240 fb _1 of data recorded by the BABAR detector at the 
PEP-II asymmetric-energy e + e~ storage rings at the Stanford Linear Accelerator Center running 
at center-of-mass energies near 10.6 GeV. 

PACS numbers: 14.40.Lb, 13. 25. Ft, 12.40.Yx 



The spectrum of known cs states can be described eco- 
nomically as two s-wave states (Df, D* + ) with J p = 
Q~, I" and four p-wave states (£>* (2317)+, L> s i(2460)+, 
£> s i(2536)+, L> s2 (2573)+) with J p = 0+, 1+, 1+, 2+, 
though the last two spin-parity assignments arc not 
firmly established. Whether this picture is correct re- 
mains controversial because the states at 2317 MeV/c 2 
and 2460 MeV/c 2 Q had been expected to lie at much 
higher masses |2|. 

We report here on a new cs state and a broad structure 
observed in the decay channels D°K + and D + K s . This 
analysis is based on a 240 fb _1 data sample recorded 
near the T(45) resonance by the BABAR detector at the 
PEP-II asymmetric-energy e + e~ storage rings. 

The BABAR detector is described in detail elsewhere . 
Charged particles are detected and their momenta mea- 
sured by a combination of a cylindrical drift chamber 
(DCH) and a silicon vertex tracker (SVT), both oper- 
ating within a 1.5 T solenoidal magnetic field. A ring- 
imaging Cherenkov detector (DIRC) is used for charged- 
particle identification. Photon energies are measured 
with a Csl electromagnetic calorimeter. We use infor- 
mation from the DIRC and energy-loss measurements in 
the SVT and DCH to identify charged kaon and pion 
candidates. 

We observe three inclusive processes Q: 



e + e 



D°K+X, D° 

D°K+X, D° 

o - 



D + K%X,D 



K~TT + 
K-TT+TT 



K° 



(1) 
(2) 

7r+7r-(3) 



For channels (1) and (2) we perform a vertex fit for the 
K~ir + and require a \ 2 probability greater than 0.1%. 
For the ir° in channel (2), we consider the photons that 
emanate from the K~ir + vertex, perform a fit with the 
7T° mass constraint, and require a \ 2 probability greater 
than 1%. The combinatorial background is reduced by 
requiring the ir° laboratory momentum to be greater than 
350 MeV/c. 

To purify the D° sample in channel (2), its quasi-two 
body decays [j| K*ir and Kp are used, allowing ranges 
of ±50 MeV/c 2 around the K* mass for Kir and ±100 
MeV/c 2 around the p mass for 7T7t. 

For channel (3), we fit two pions with the same charge 
and a kaon of the opposite charge to a common vertex 
to form the D + candidate, and require a \ 2 probability 
greater than 0.1%. We obtain the Kg sample with a 
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FIG. 1: (a) K~n + , (b) R-tt+tt and (c) K~-k + -k + mass 
distributions for all candidate events to channels (1), (2) and 
(3) respectively. The shaded regions indicate the definition of 
signal and sidebands regions. 



fit that constrains the mass and require a x 2 probability 
greater than 2%. Kg candidates are retained only if their 
decay lengths are greater than 0.5 cm. 

For all three channels, the D candidate is combined 
with an identified K, requiring a vertex fit \ 2 probability 
greater than 0.1%, and constraining the vertex to be in 



the i 



luminous region. To reduce combinatorial back- 



ground from the continuum (e H 



qq,q = u,d,s,c) 



and B-meson decays, each DK candidate must have a 
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FIG. 2: The DK invariant mass distributions for (a) D° K _ 7T+ K+, (b) D° K _ n+iT oK+ and (c) D+_^ + ^ + K°. The shaded 
histograms are for the D-mass sideband regions. The dotted histogram in (a) is from e + e~ — ^ cc Monte Carlo simulations 
incorporating previously known D s states with an arbitrary normalization. The insets show an expanded view of the 2.86 
GeV/c 2 region. The solid curves are the fitted background threshold functions from the three separate fits described in the 
text. 



momentum p* in the e + e center-of-mass frame greater 
than 3.5 GeV/c. 

Figs. Ufa), (b) and (c) show the R-tt+, K-jt+tt , 
and K~tt + tt + invariant mass distributions, respectively. 
All distributions show pronounced peaks at the D mass, 
with signal yields of about 950,000, 790,000, and 430,000 
events. Fits using a polynomial and a single Gaussian 
give a = 7.6,12.6,6.0 MeV/c 2 for the three widths. We 
define the signal region by ±2er from the fitted D mass 
and establish sidebands at (— 6er, — 4cr) and (4er, 6er). In 
the signal regions, the signal-to-background ratios are 
4.1, 1.2, and 2.2 respectively. 

Selecting events in the D signal regions, Fig. [3 shows 
the D°K + invariant mass distributions for channels (1) 
and (2), and the D + Kg invariant mass distribution for 
channel (3). To improve mass resolution, the nominal 
D mass and the reconstructed 3-momentum are used to 
calculate the D energy for channels (1) and (3). Since 
channel (2) has a poorer D° resolution, each K~ir + -K a 
candidate is kinematically fit with a D° mass constraint 
and we require a % 2 probability greater than 0.1%. 

We find that the fraction of events having more than 
one DK combination per event is 0.9% for channels (1) 
and (3) and 3.4% for channel (2). In the rest of the paper, 
we use the term reflection to describe enhancements pro- 
duced by two or three body decays of narrow resonances 
where one of the decay products is missed. 

The three mass spectra in Fig. [21 present similar fea- 
tures. 



• A single bin peak at 2.4 GeV/c 2 due to a reflection 
from the decays of the D sl {2536) + to D*°K+ or 
D* + Kg in which the ir° or 7 from the D* decay is 
missed. This state, if J p = 1 + , cannot decay to 
DK. 

• A prominent narrow signal due to the _D S 2(2573) + . 

• A broad structure peaking at a mass of approxi- 
mately 2.7 GeV/c 2 . 

• An enhancement around 2.86 GeV/c 2 . This can 
be seen better in the expanded views shown in the 
insets of Fig. 

In the following we examine different background 
sources: combinatorial, possible reflections from D* de- 
cays, and particle misidentification. 

Backgrounds come both from events in which the can- 
didate D meson is correctly identified and from events 
in which it is not. The first case can be studied combin- 
ing a reconstructed D meson with a kaon from another 
D meson in the same event, using data with fully recon- 
structed DD pairs or Monte Carlo simulations. No signal 
near 2.7 or 2.86 GeV/c 2 is seen in the DK mass plots for 
these events. The second case can be studied using the D 
mass sidebands. The shaded regions in Fig. |3 show the 
DK mass spectra for events in the D sideband regions 
normalized to the estimated background in the signal re- 
gion. No prominent structure is visible in the sideband 
mass spectra. 
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We examined the possibility that the features at 2.7 
and 2.86 GeV/c 2 could be a reflection from D* or other 
higher mass resonances. Candidate DK pairs where the 
D is a D*-decay product are identified by forming Dir 
and D7 combinations and requiring the invariant-mass 
difference between one of those combinations and the D 
to be within ±2cr of the known D* — D mass difference. 
No signal near 2.7 or 2.86 GeV/c 2 is seen in the DK mass 
plots for these events. Events belonging to these possi- 
ble reflections (except for the D*° — > D 7 events, which 
could not be isolated cleanly) have been removed from 
the mass distributions shown in Fig. (corresponding to 
«8% of the final sample). 

We use a Monte Carlo simulation to investigate the 
possibility that the 2.7 or 2.86 GeV/c 2 signals could be 
due to reflections from other charmed states. This simu- 
lation includes e + e~ — > cc events and all known charmed 
states and decays. The Monte Carlo events were gener- 
ated using a detailed detector simulation and subjected 
to the same reconstruction and event-selection procedure 
as was used for the data. The D°K + effective mass dis- 
tribution for these Monte Carlo events is shown (dotted) 
in Fig. |2fa) for channel (1). The normalization is ar- 
bitrary. No peak is found in the 2.7 and 2.86 GeV/c 2 
D°K + signal regions. We note that the simulation un- 
derestimates the size of the Z? s i(2536) + reflection and the 
D S 2(2573) + signal relative to the background. No such 
discrepancy is found in the study of the D°tt + final state, 
therefore we attribute this effect to a poor knowledge of 
the strange-charmed meson cross sections. In order to 
improve the data-Monte Carlo comparison, the events 
having a Z? s i(2536) + and _D S 2(2573) + in the final state 
have been scaled up by factors 5 and 2 respectively. 

We checked the possibility that the structures at 
2.7 and 2.86 GeV/c 2 are due to misidentifying pions 
as kaons by assigning the kaon mass to the pion in 
D°tt + data events. We observe no structure near 2.7 or 
2.86 GeV/c 2 in the resulting D°K + invariant mass dis- 
tribution. Monte Carlo simulations and tests using the 
data show that these structures also do not originate from 
protons misidentified as kaons from high mass charmed 
baryon decays. 

Wrong sign D°K~ mass distributions for channels (1) 
and (2) have also been examined and we find no signal 
in either mass spectrum. 

A more detailed study in channel (1) of the 2.7 GeV/c 2 
structure shows a broad structure in this mass region for 
events from the D° sidebands in which the DK candidate 
has a very low p* (p* < 3 GeV/c). This is not seen 
in channels (2) or (3) however. We conclude that the 
assignment of the 2.7 GeV/c 2 structure to a reflection 
remains inconclusive. 

By comparing the reconstructed mass distributions for 
the DK system with those generated with Monte Carlo 
simulations, we obtain the mass resolutions. The resolu- 
tions are similar in the three channels, increasing linearly 



from 1.7 MeV/c 2 at a mass of 2.5 GeV/c 2 , to 3.5 MeV/c 2 
at a mass of 2.86 GeV/c 2 . 

In the following discussion we label as Z? s ,/(2860) + 
the structure in the 2.86 GeV/c 2 mass region and as 
X(2690) + the structure observed in the 2.7 GeV/c 2 mass 
region. We fit to the three DK mass spectra shown 
in Fig. H from 2.42 GeV/c 2 to 3.1 GeV/c 2 (exclud- 
ing the Z3 s i(2536) + reflection) using a binned \ 2 min- 
imization. The background for the three DK mass dis- 
tributions is described by a threshold function: (m — 
m t h) a e -P m -~< m2 - 5m3 where m th = m D + m K . A fit 
to the Monte Carlo distribution shown in Fig. |2Ia) us- 
ing this background expression and one spin-2 relativis- 
tic Breit-Wigner for the £> S 2(2573) + gives a good fit 
with 32 % x 2 probability. In the fit to the data, the 
D S 2(2573) + and D s j(2860) + peaks are described with 
relativistic Breit-Wigner lineshapes where spin-2 is as- 
sumed for the .D S 2(2573) + and spin is used for the 
L> sJ (2860)+. We find that the L> sJ (2860) + parameters 
are insensitive to the choice of the spin. The best de- 
scription of the X(2690) + structure is obtained using a 
Gaussian distribution. The results from the fits are sum- 
marized in Table [Q Table |n| summarizes the \ 2 proba- 
bilities, the number of D s j (2860) + events (with statisti- 
cal and systematic errors) and the D s j(2860) + statistical 
significances from the three separate fits to the DK mass 
spectra. 

The fits give consistent values for the parameters of 
the three structures. We notice a smaller width of the 
D s2 (2573) + in the D^_^ +7T+ K® channel which we at- 
tribute to the uncertainty in the description of the back- 
ground. We compute also the ratios of the yields of 
,D s j(2860) + with respect to Z? s2 (2573) + finding agree- 
ment, within statistical errors, between the three chan- 
nels. 

The presence of resonant structures can be visually 
enhanced by subtracting the fitted background threshold 
function from the data. Fig. 01 shows the background- 
subtracted D° _ + K+, D%_ , K+, and D+ , ,K° 
invariant mass distributions in the 2.86 GeV/c 2 mass re- 
gion. Fig. EJd) shows the sum of the three mass spectra. 

We also fit to the three distributions simultaneously. 
The parameters from this fit are labelled DK A in Table[Q 
If we remove the £> s j(2860) + , the \ 2 increases by 108 
units while the number of degrees of freedom increases 
by five. 

As a systematic test, we repeated the fits varying the 
lower p* cut on the DK system from 3.50 to 3.75 and to 
4.00 GeV/c. We also restricted the fit to the D sJ (2860)+ 
only and replaced the threshold function which represents 
the background with a polynomial. Fits have also been 
performed without removing the events associated to D* 
reflections and modifying the spin of D S 2(2573) + . The 
systematic uncertainties take into account the variation 
of the resonance parameters among the three different 
final states and the resonance parameterizations. The 
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m(D K) GeV/c 2 



FIG. 3: Fitted background-subtracted DK invariant mass distributions for (a) D° K ^^ + K + , (b) D ( ^ c ^ 7r+7l0 K + , (c) D~^_ + + K%, 
and (d) the sum of all modes in the 2.86 GeV/c 2 mass region. The curves are the fitted functions described in the text. 



TABLE I: Results from the fits to the total DK mass spectra of Fig. [5] Quantities are in units of MeV/c 2 . Errors are statistical 
only. Simultaneous fits of the three mass spectra are labelled with DKa and DA'a. 

Fit m(P s2 (2573)+) r(P s2 (2573) + ) m(X(2690)+) cr(A(2690)+) m(P sJ (2860) + ) r(P sJ (2860) + ) X 2 /NDF 

D u +K + 2572.4 ± 0.4 27.6 ± 0.5 2687 ± 4 41 ± 5 2855.4 ± 2.0 37 ± 8 26/20 

D K-l+n» K+ 2572.3 ± 0.5 28.4 ± 0.9 2682 ± 5 52 ± 5 2860.8 ± 4.0 52 ± 14 33/20 

D K-*+-x+ K ° 2572.6 ± 0.9 21.9 ± 1.1 2684 ± 7 50 ± 7 2856.6 ± 8.0 81 ± 25 26/21 
DKa 2572.3 ± 0.3 27.1 ± 0.6 2684 ± 3 48 ± 2 2856.6 ± 1.5 47 ± 7 100/72 

Fit m(P 32 (2573)+) r(P s2 (2573) + ) m(X(2690)+) r(A(2690)+) m(P sJ (2860) + ) r(P sJ (2860) + ) X 2 /NDF 

DK B 2572.3 ± 0.3 27.0 ± 0.5 2688 ±4 112 ± 7 2857.6 ± 1.9 38 ± 7 112/72 



TABLE II: x 2 probabilities, P s j(2860) + event yields and sta- 
tistical significances from the three separate fits to the total 
DK mass spectra of Fig. 



Channel 


x 2 


D sJ (2860) + 


Statistical 




probability (%) 


events 


significance 


n u k+ 

U K-TT+ K 


17 


886 ± 134 ± 49 


6.2 cr 


D° K -, + ,oK+ 


3 


1146 ± 157 ± 78 


6.5 cr 




21 


371 ± 84 ± 53 


3.7 a 


DK A 


1.6 


2717 ± 262 ± 190 


8.4 cr 


DK b 


0.2 


2161 ± 238 ± 151 


7.7 cr 



uncertainty on the mass scale is estimated to be of the 
order of 1 MeV/c 2 . 

We obtain the mass and width of D S 2(2573) + : 

m(D s2 (2573)+) = (2572.2 ± 0.3 ± 1.0) MeV/c 2 
r(D s2 (2573) + ) = (27.1 ± 0.6 ± 5.6) MeV/c 2 , 

where the first errors are statistical and the second sys- 



tematic. For the new state we find 

m{D sJ (2860)+) = (2856.6 ± 1.5 ± 5.0) MeV/c 2 
r(L> sJ (2860) + ) = (47 ± 7 ± 10) MeV/c 2 . 

Since the assignment of the X(2690) + as a reflection 
is inconclusive, the three mass spectra have also been 
fit including the Jf(2690) + as an additional resonance 
(Breit-Wigner, rather than Gaussian shape). This gives 
fit DKb shown in Table HJ The resulting resonance pa- 
rameters are: 

m(X(2690)+) = (2688 ± 4 ± 3) MeV/c 2 
T(X(2690) + ) = (112±7±36) MeV/c 2 . 

In summary, in 240 fb _1 of data collected by the 
BABAR experiment, we observe a new D+ state in the 
inclusive DK mass distribution near 2.86 GeV/c 2 in 
three independent channels. The decay to two pseu- 
doscalar mesons implies a natural spin-parity for this 
state: J p = + , 1~, . . . . It has been suggested that this 
new state could be a radial excitation of Z?*j (2317) 
although other possibilities cannot be ruled out. In the 
same mass distributions we also observe a broad enhance- 
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ment around 2.69 GeV/c 2 which it is not possible to as- 
sociate to any known reflection or background. 
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